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ABSTRACT 


The purpose of this investigation was to try to develop 
a technique for the remote sensing of water resources by in- 
ferring river discharges from power spectral analyses of river 
meander patterns. We used ERTS-1 data to help determine whether 
a correlation between river meander power spectra and river 
discharge frequency exists. In the course of this study we 
developed techniques for analysing remote sensing imagery of 
river channel patterns; we discovered an important new re- 
lationship between the discharge spectra of rivers and the 
time dependence of their flood recessions, which is important 
to hydrologic modeling of stream flow; and, within the un- 
certainty imposed by the data we were able to analyze, we 
showed that only a rough correlation exists between the meander 
power spectra and the discharge spectra of rivers, and that 
other variables must also influence this relationship pre- 
venting a simple inversion of the meander spectrum to obtain 
a discharge spectrum. 
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power law depender. ,o of the discharge frequency 
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Figure 10. Power law fits to daily discharge data (circles) 
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1. Introduction 

1.1 Background 

A large number of correlations between some average 
meander wavelength and some characteristic discharge have 
been proposed [see for example Jefferson, 1902? Inglis, 1949; 
Leopold and Wolman, 1957; Dury, 1965? Carlston, 1965; Schumm, 
1971]. These correlations clearly suggested that there was 

t 

a relationship between meander wavelength and discharge but 
they failed to agree on its quantitative form. Some of this 
disagreement results from the oversimplification inherent in 
using a single meander wavelength and a single discharge to 
characterize the river rather than using the complete spectra 
of wavelengths and discharges. 

Speight (1965, 1967) appreciated that the entire oscilla- 
tory pattern of a river must be iirportant in characterizing 
its meandering and presented power spectra of the auto-correla- 
tions of the directions of flow measured at equally spaced 
points on the talwegs of several Australasian rivers. These 
meander power spectra showed structure which we interpreted 
as an indication that several characteristic length scales may 
be required to quantitatively describe a meander pattern. The 
idea of using a spectral analysis of the reach of a river as 
the basis of a correlation rather than a subjective estimate 
•.,f an assumed single length scale is a necessary generaliza- 
tion in describing the connection between a rivet's meander 


pattern and its discharge. However/ he retained the idea 
that a single discharge could be correlated with the multi- 
ple length scales. Just as there is an essential difficulty- 
in attempting to characterize a meander pattern by a single 
length scale, there is a fundamental problem in trying to 
choose the dominant discharge, i.e. that discharge most 
effective in establishing the system of meanders. We in- 
vestigated the possibility that the further generalization 
of the correlation to include the time-behavior of the dis- 
charge might bring an order to the' relationship between the 
total meander pattern and the complete record of the discharge. 
This more general correlation should be sufficiently reliable 
to quantitatively assess a river's flowrate from a spectrum 
of its meanders, thus making the knowledge of a region's 
water resources accessible from aerial or satellite imagery 
of the area. 

1.2 Outline of Study 

Our study of a possible correlation between the 
stream meander power spectrum and the stream discharge fre- 
quency distribution required the accomplishment of the 

m 

following tasks: 

1 . Selection of appropriate rivers . 

2. Collection of suitable photographic coverage. 

3. Collection of historical streamflow data. 

4. Digitization of streamflow data. 

5. Digitization of stream meander patterns. 
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6. Matching of ' individually digitized portions of 
meander patterns to obtain a continuous record - 

7. Construction of stream discharge frequency dis- 
tribution from discharge data. 

8. Construction of power spectra of local river 
meander directions. 

9. Determination of relationship between meander 
power spectra and discharge frequency distributions. 

The completion of the first eight tasks for the Feather 
River in California and a number of other rivers throughout 
the United States, established general procedures for generat- 
ing the stream discharge probability density functions and 
meander power spectra for all rivers. 

From these analyses we have shown that a rough relation- 
ship appears to exist between the discharge spectra and meander 
spectra of most of the rivers studied. From the data for this 
study we were not able, however , to determine a qualitative 
relationship suitable for uniquely determining the discharge 
spectrum of a river from remote sensing imagery of its meander 
pattern. But such a relationship may yet be realizable from 
a broader data base. * 

1.3 Data Acquisition 

Prior to attempting to obtain the imagery and 
hydrologic data required for our study, we first had to 
develop selection criteria for defining rivers suitable 


for study, 
1 . 


These were 


Uniformity of geology . 

2. Length of reach with minimum spatial variation in 
discharge axong the course. 

3. Length of reach with minimum seasonal variation in 
discharge . 

4. Resolvability of river on ERTS-1 imagery and/or 
availability of historical photographic coverage. 

5. Availability of continuous historical hydrologic data 

Using these criteria we selected sixteen rivers which 

sampled a wide range of discharge rates and a broad geographical 
distribution. 


TABLE 1 


LIST OP RIVERS STUDIED 

Median 


River 

Imagery 

Hydrologic Station 

Q 

ft. 3/sec 

Length of Reach 
(ft.) 

Mississippi River 

ERTS-1 

7/0320,2890 

350,000 

4.1xl0 6 

Sacramento River 

ERTS-1 

11/3780,3890,3895 

8,000 

6, 0x10 5 

Feather River 

ERTS-1 & Aerial 

11/4070 

3,000 

1.2xl0 5 

Manistee River 

Aerial 

4/1235,1240,1260 

200 - 1000 

5. 0x10 5 

Chippewa River 

Aerial 

5/3565 

1,000 

1.8x10 5 

Ontonagon River 

Aerial 

4/0400 

900 

7 . 5x10 5 i 

Flambeau River 

Aerial 

5/3585 

800 

1.3xl0 b . 1 

Mushegon River 

Aerial 

4/1215 

700 

1.2xl0 S 

Manistigue River 

Aerial 

4/0495,0550 

400 - 600 

1.6x1 0 5 

Homochit to River 

Aerial 

7/2910, 2925, 2945 

400 

2.6xl0 5 

Red River 

Aerial 

5/0540 

300 

3.9xl0 5 

Bad River 

Aerial 

4/0270 

250 

1.3xl0 S 

White River 

Aerial 

4/0275 

200 

l.lxlO 5 

Boise Brule River 

Aerial 

4/0255 

150 

1.3xl0 5 

Rifle River 

Aerial 

4/ 1395,1405,1420 

80 - 200 

1. 7xl0 5 

Jump River 

Aerial 

5/3615,3620 

70 - 140 

1.8xl0 5 
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1.3.1 ERTS-1 Imagery and Aerial Photography 

Of these sixteen rivers (see Tablt> 1) only three- 
were of sufficient size that their meander patterns were 
clearly resolvable in the IR imagery from ERTS-1. These 
were, the Feather and Sacramento Rivers in California, and the 
Mississippi River from Tennessee to New Orleans. These rivers 
also had largest mean discharges (see Table 1} . Fourteen 
other rivers, which were needed to broaden the range of dis- 
charges but were not resolva ble in the ERTS-1 imagery, had 
to be studied using aerial photography. Table 1 lists all of 
these rivers giving the type of imagery used, the hydrologic 
stations from which discharge data was available, the median 
discharge and the total length of the reach of river studied. 
Where a range of discharges is listed, corresponding to 
values at different stations, shorter reaches associated with 
each station were analyzed. We have obtained aerial photo- 
graphs (both infrared and panchromatic) of rivers from the 
Agricultural Stabilization and Conservation Service of the 
U.S. Department of Agriculture, Department of Water Resources 
of the California State Resources Agency, the Topographical 
Division of the U.S. Geological Survey, and the Cartographic 
and Audiovisual Records Division of the National Archives and 
Records Service. Infrared satellite imagery from ERTS-1 has 
been provided by the NASA Data Processing Facility at Goddard 
Space Flight Center. 
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1.3.2 Hydrologic Data 

The streamflow data for all of the rivers used 
in this study consisted of daily discharge measurements 
from gauging stations operated by the TJ.S- Geologic Survey'. 
This data was obtained in machine measurable format from the 
Water Resources Division of the U.S. Department of the 
interior. The discharge records for some of these rivers, 
such as the Sacramento, encompass as much as seventy five 
years of continuous measurement. 

2. Accomplishments 

As a basis for this study we have developed a fully auto- 
mated system for obtaining both the discharge and meander 
wavelength spectra. Discharge spectra (probability of dis- 
charge per unit discharge vs. discharge) are constructed from 
historical records of daily stream discharge. Generation of 
meander power spectra involves three elements: digitization 

by photoelectric optical tracking of stream banks on each 
frame of photographic or television imagery? collation and 
matching of successive frames into a single data record for 
each stream; and a Fourier transform analysis of the data. 

This system has been developed to facilitate the analysis of 
the large number of rivers required to assure the statistical 
reliability of the correlation. 

From the imagery of a selected river reach, we determine 



the positions of the curve with respect to a cartesian coordinate 
system. This description of the river's course is transformed 
to a (9,s) representation, where 0 measures the angle that 
the curve makes with a reference direction as a function of 
distance travelled along the curve s. The (9,s) description 
is preferred since the meander pattern of a river may be 
represented by a multivalued function of position in a cartesian 
representation. The meander power spectrum which we calculate 
is the power spectral density of 9(s). 

2.1 Digitization 

The digitization of river meander patterns from 
aerial imagery is most accurately and economically accomplished 
through photoelectric optical scanning. We have developed a 
program for digitizing river meanders, using commercially 
available machines employing this technique. An important 
condition on the digitization procedure is that data points 
be located at equal increments of distance along the meander 
curve. This condition follows from the fact that local meander 
direction, 9, is a function of distance along the meander and 
the algorithm used for constructing the power spectrum 
requires that we know this quantity at equal increments of 
distance » 

The essence of the digitization procedure is as follows 
(see also Fig. 1) : The initial point on the meander pattern 

is found by scanning along a horizontal or vertical line and 



measuring the optical density profile along the scan (line AG 
in Pig. 1). The river bank, i.e., the point digitized, is 
defined as the location of the point of maximum gradient in 
the density profile (point 1) . The second point is determined 
by an iterative process starting with a scan (line CD) parallel 
to the initial scan but displaced by a distance, s, from point 
1. The first estimate of this point (2 # ) is determined in the 
same manner as above. The distance between points 1 and 2 ' 
is then calculated and if it is not equal to the required 
spacing, s, plus or minus some small A, another scan (line EF) 
is made alone a line perpendicular to that line connecting 
points 1 and 2 ; at a distance s from point 1. Point 2 ° is 
then determined along this scan. If the distance between it 
and point 1 is still not within s ± A, the iterative process 
is repeated until convergence is obtained. Once point 2 is 
located, the search for point 3 begins along a scan line (GH) 
perpendicular to the line connecting points 1 and 2 at a distance 
s from point 2. In this manner the machine proceeds along the 
meander curve digitizing points at equal distance increments 
along the curve. 

As is generally ,the case, the imagery of the river con- 
sists of a number of overlapping frames, thus this digitization 
procedure is repeated for each frame. The data for adjacent 
frames must then be matched to give a continuous digitized 
record of the meander pattern. Because of the large overlap 
between frames (roughly one-third of the data on the end of 
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each frame overlaps with the data on the beginning of the 
next frame) the data sets can be uniquely matched for con- 
gruency in the overlap region. We have developed a computer 
program which finds the appropriate coordinate transformation, 
i.e., includes both translation and rotation of one frame with 
respect to the other. This is accomplished by considering a 
length of river roughly half the length of the overlap region 
on one frame and effectively sliding this portion of the data 
along the overlap portion of the adjacent frame, finding that 
transformation within the overlap region which minimizes the 
sum of the squares of the distances between matched points. 

Once the appropriate coordinate transformation is determined 
rail of the data points on the second frame are transformed. 

The process is then repeated to match successive frames 

until the entire record is transformed into a single coordinate 

system. 

For purposes of comparison a computer generated plot of the 
digitized meander pattern of a section of the Feather River 
is shown as a slightly displaced overlay in Figure 2 together 
with the infrared photograph from which the pattern was derived. 
As can be seen this technique provides an accurate reproduction 
of the meander pattern. A sample of the final product of this 
digitization and matching procedure is shown in Figure 3 for 
a 23 mile reach of the Feather River below the dredge tailings 
at Oroville to the Southern Pacific Railroad bridge north of 



- ■ - 
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Figure 2. Tnf rargd imagery of section of Feather River with Slightly displaced 
overlay of digitized meander pattern 
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Yuba City. A topographic map of this reach of river is shown 
in Figure 4. 

2.2 River Meander Power Spectra 
2.2.1 Technique 

The digitization and matching procedures 
described above produce a set of data points (S,Y coordinates) 
which are equally spaced along the course of the river. A 
power spectral analysis of the river cannot be made directly 
from the X versus Y data since the river may double back upon 
itself making X a double valued function of Y. An equivalent 
representation of the river, which is single valued and thus 
amenable to power spectral analysis is its local direction, 

0, as a function of the distance, s, along the river's course. 

2 3 1 * ■ * 

The power spectral density (deg /10 ft.) for the direction 
0 is computed using standard techniques for determining the 
autocorrelation function, smoothing, and taking the Fourier 
transform (e.g., J.S. Bend at and A.G. Piersol, Measurement and 
Analysis of Random Data , John Wiley, 1966) . 

These spectra are calculated as a function of wave 
3 

number (per 10 ft. ) . The range of wave number is limited at 
the high end (short wavelength) by the interval between data 
points and at the low end (long wavelength) by the length of 
the reach Fourier analyzed and the number of degrees of 
freedom. The maximum wave number is given by the Nyquist 
criterion to be N/2-t where N is the number of data points 


and £ is the length of the reach. The minimum wave number 
at which spectral information is obtained is N/2-C.m where m 
is the number of spectral estimates, or alternatively n/4-t 
where n — 2N/m is the number of degrees of freedom, a measure 
of the uncertainty in the power spectral density. The 
analyses are all for 6 or more degrees of freedom which gives 
us 80% or more confidence that the spectral value lies within 
2.8 and 0.5 times the computed estimate. 

2.2.2 Characteristics of Spectra 

The meander spectra which we have calculated 

for the various rivers described above are all rather simple, 

broad spectra having significant power over a wave number 

range of more than a decade. All of these spectra are shown 

in Appendix I. Three examples are also shown in Figuure 5, 

together with scaled traces of the meander patterns which have 
% 

been spectrally analyzed. The error bar on each spectrum 
indicates the 80% confidence interval for the power spectral 
density determined from the number of degrees of freedom 
(Blackman and Tukey, 1958). 

At intermediate wave numbers each of these spectra have 
an apparent power-law dependence of the spectral density on 
wave number, as shown by the dashed line preliminary fit to 
the spectral estimates. The exponent of the power law is 
given by the indicated slope which is probably accurate to 
within ± 20%. At lower wave numbers the spectrum flattens. 


suggesting a broad peak at a power density two to three 

orders of magnitude above the noise level. The wave numbers 

at which the power law dependences break to flatter dependences 

on the Ontonagon, Homochit to and White Rivers, correspond to 

3 4 3 

wave lengths of 5x10 , 10 and 2.5x10 ft., respectively and 
reflect recognizable scales in the meander patterns of each 
river shown in the figure. At high wave numbers the power 
law dependence is truncated at the relatively constant level 
of the noise which extends up to the Nyguist limit on the 
wave number- Though there are small scale fluctuations about 
the general spectral structure outlined above, they are not 
significant within even the 80% confidence interval. 

Some streams, however, do show more complex structure as 
can be seen in Fig. 6. Clearly the off-set structure shown 
’ in the spectrum of the Manistee River is significant within 
the 80% confidence interval. There is also a suggestion of 
similar structure in the spectra of the Bad and Bois Brule 
Rivers though not at the same level of confidence. 

Thus the spectra which we have investigated so far can 
be described by one or more linear power law segments of 
differing slopes and magnitude which break at characteristic 
wave numbers. 

We have investigated the stationarity of meander power 
spectra by constructing spectra for a series of subreaches of 
the Manistee River to study the variation of its spectrum 
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Figure 5. Pcwer spectral density (deg 10 ft.) ‘versus wave number 
(per lO^ft.) and wave length (IQ^ft.) for three river readies shaving 
rather simple, broad spectra with power law dependent segments 
approximately fitted by dashed lines with slopes indicated. A scaled 
•trace of the meander pattern of each reach analyzed is also shown, 
together with data on the length of the reach, the number of data 
points., . and the nurrber of spectral estimates, and an error bar 

indicating the 80% confidence internal for the power spectral density 
estimate. 
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MEANDER POWER SPECTRAL 



Figure 7. Stationarity of the offset structure in the power spectrum of the 
Manistee River shown by its persistence in the spectra of four separate 
readers spanning a total length of sixty miles. 









along the course of the stream. Figure 7 shows spectra for 

4 

four such consecutive subreaches, each of length 7.3x10 ft. 
at the downstream end of the larger reach of the Manistee River 
shown in Figure 6. Each subreach may be recognized in the 
larger reach of the Manistee River by the traces of their 
meander patterns . As can be seen in Figure 7 the spectra of 
the four subreaches are quite similar; the slopes and magni- 
tudes of the power- law segments of the spectra and the wave 

-1 -1 

numbers of breaks in the spectra at 10 and 5x10 are the 
same in all four spectra to well within the 80% confidence 
interval. Thus these characteristic features of the power 
spectrum of the Manistee River are stationary over a length 
of at least 60 miles along its course. We have also con- 

4 

structed the power spectrum of a 7.3x10 ft. subreach of the 
Manistee River located beyond the upstream end of the reach 
shown in Figure 6. This spectrum has only a single power 
law dependent segment, breaking at a wave number of 1 and 
differing significantly from those shown in Figure 7. This 
difference is not surprising, however, in view of the fact 
that the median discharge in this upper subreach also differs 
by a factor of 5 from that in the lower subreaches. 

In comparing our meander power spectra with those of 
earlier investigators (Speight, 1965, 1967? Toebes and Chang, 
1967) we note first that they used linear scales to plot the 
power spectral density estimates- Thus they did not observe 
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the power law segments of meander spectra identified in this 
paper. Speight (1965, 1967) appears to have placed an unwarranted 
credence in the physical significance of a great many “spectral 
peaks", nearly all of which we would classify as random fluctua- 
tions on the power law portions of the spectra. In fact none 
of Speight's "peaks" were resolvable on the basis of the 
number of degrees of freedom which he used in his spectral 
analysis, but he attributed significance to them on the basis 
of some stationarity over the reach. On the other hand, 

Toebes and Change (1967) have gone to the opposite extreme in 
suggesting that meander spectra are inherently nonstationary, 
reflecting only a randomness in the meander patterns. Neither 
of these investigators analyzed a large enough number of 
meander power spectra to justify their respective views. 

Based on the rivers which we have studied so far (examples 
of which have been presented here) we believe that there is 
significant structure in a meander power spectrum, namely the 
slopes and magnitudes of the power law segments and the wave 
numbers at which breaks in these segments occur. It is these 
characteristic parameters of the meander spectra which we will 
attempt to correlate with such characteristics of the discharge 
spectra as the modal discharge and the exponent of the flood 
recession. 

The cause of this simple power law relationship is not 
yet understood but because of the generality of the relationship 





it may be fundamental for the process of meander formation. 

2.3 Discharge Spectra 

2.3.1 Technique 

The discharge spectrum or frequency distribu- 
tion of a river is the probability per unit discharge that 
its discharge, or flow rate, Q lies within the interval Aq at 
Q. The integral of this distribution over discharge is the 
flow duration curve commonly used in hydrologic studies (see, 
for example, Chow, 1964) . From an historical discharge 

record of length T, the frequency distribution which we calcu- 
1 At 

late is ~ . where At is the portion of the time the dis- 

charge lies within Aq at Q. 

2.3.2 Power Law Dependence 

From historical records of daily river dis- 
charge, we have constructed the frequency distributions from 
the fraction of time the discharge lies within a prescribed 
interval per unit interval. These are all shown in Appendix I. 
A remarkable property of most of these discharge frequency 
distributions is their nearly linear character on log-log plots 
for values of discharge larger than the mode (Schubert and 
Lingenfelter , 1974) . Several examples of this are shown in 
Figure 8. The daily discharge data on which each distribution 
was based extends over the indicated time interval for the 
particular gauging station identified by number according to 
the convention adopted by the U.S.G.S. (1964) . The slopes 
of the linear portions of these distributions vary from river 
to river and even from station to station on the same river 
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Red River 


Bad River 


STATION 5/3620 1915-1970 STATION 7/3370 . 1936-1970 STATION 4/0270 1948-1970 













over the range from less than -1 to -5. Plow duration 
curves have previously been interpreted as representing a 
random process described, for example, by a log-normal dis- 
tribution (Chow, 1964; Leopold et al., 1964). Our results, 
however, show that such an interpretation is not appropriate 
for a large number of rivers, since it is inconsistent with 
the clearly linear character of the log- log frequency dis- 
tribution plots. Out interpretation of this relationship 
is discussed in section 2-4. 

2.4 Relationship Between Discharge Spectra and River 
Flood Decay 

We have shown (Schubert and Lingenfelter, 1974) 
that the distribution for discharges greater than the mode 
is not stochastic as previously suggested but is essentially 
deterministic in nature, reflecting the decay phase of the 
flood hydrograph. We have not found any previous suggestion 
of a direct relationship between the form of the flow dura- 
tion curve and the flood hydrograph. 

If s is the slope of the log-log frequency distribution, 

then 

^ * Q S - (2.4.1) 


Integrating eq. (1) we find 
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1 


Q « t 


s+1 


( (2.4.2) 


Since the observed values of s lie between about -1.5 and -5 
the exponent would range from about -2.0 to -0.25. We 

suggest that equation (2.4.2) represents the decay phase of 
the flood hydrograph, where t is measured from a time tg near 
the flood peak. The time tg can be uniquely determined from 
any two discharge measurements Q^, Q 2 at times t^, t 2 during 
the flood recession by 



The inverse power law dependence of the discharge on time, 
which we find here, differs from the superposition of several 
exponential decay curves, which have previously been used 
(Chow, 1964) to empirically fit the flood recession. 

To test our suggestion that the linear nature of many 
of the log-log frequency distributions represents the recession 
portion of the flood hydrograph, we have compared the time 
dependence of the discharge predicted by equations (2.4.2) 
and (2.4.3) with the measured decay of discharge following 
individual flood peaks on the various rivers studied. We find 
that the predicted decay at each station does indeed describe 
the measured discharge following all flood peaks at that 
station. A typical example of 4 :h« agreement between predicted 


and measured flood recession is shown in Figure 9 for the 
Sacramento River near Red Bluff, California, in 1936. As 
can be seen, the curves of the theoretical flood decay are an 
excellent fit to the data points which indicate the measured 
values of the daily discharge. The theoretical curves are 
based on a value of s equal to -2.6 (see Figure 8} and 
values of tg equal to 16 January, 22 February and 4 April 
for the respective floods shown in Figure 10. At this 
station on the Sacramento River floods decay according to 
the rule t ^.625^ allows the recession to be determined 

for as long as a month following the flood peak. From the 
hydrograph of Figure 10 it can be seen that the deterministic 
flood decay extends down to discharges of about 10 c.f.s. at 
this station. Belov/ this discharge level the flow rate 
variations appear to be stochastic in nature. This also is 
consistent with the fact that the linear relationship in the 
log -log frequency distribution (Figure 9) ceases at dis- 
charges below about 10^ c.f.s. at this station. 

Additional examples of the agreement between predicted 
power law flood decays and measured flood recessions are shown 
in Figure 10 and 11 for 3 consecutive flood decays on the 
Homochitto River (7/2910) and the Bad River (4/0270) , re- 
spectively. The power law decays are excellent fits to the 
observed values of daily discharge during the flcod recessions. 
Values of s equal to -1.9 and -1.8 (from Figure 8) for the 



Homochitto and Bad Rivers, respectively, were used to con- 
struct the power law decays shown in Figures 10 a- i 11. The 

floods at this station on the Homochitto River decay according 
-1 11 

to t , while on the Bad River station the recession 

follows the law t 

It is noteworthy that the flood decays on the Sacramento 

River (F\gure 9) extend from discharges of about 10^ c.f.s. 

down to 10 4 c.f.s., while on the Bad River (Figure 11) the 

recessions extend from 10 c.f.s. down to below 10 c.f.s. and 

on the Homochitto River the transients cover the discharge 

3 2 

range from more than 10 c.f.s. to below 10 c.f.s. The power 

law flood decays are thus seen to be excellent representations 

of the measured discharges for recessions which extend over 

3 orders of magnitude in discharge.- Moreover these rivers are 

from widely separated regions and represent widely differing 

hydrologic systems. From Figure 8 it can be seen that the 

Red River (5/0540) is similar to the Sacramento in that power 

5 4 

law flood decays range from 10 c.f.s. down to below 10 c.f.s. 
Also from Figure 8 we note that the Manistee River (4/1260) 
like the Baci River (4/0270) has power law flood decays in the 
discharge range 10 4 - 10^ c.f.s. and the Pine (4/1250), Thunder 
Bay (4/1325) and West Branch Manistigue (4/0550) Rivers, like 
the Homochitto River (7/2910) , have power law flood recessions 
in the discharge range 10 - 10 c.f.s. The Jump River (5/3620) 

has power law transients over a two order of magnitude discharge 




Jan. Feb. Mar. Apr. 

1936 

Figure 9. Typical daily discharge data (circles) showing the recessions of 
three floods on the Sacramento River near Red Bluff in 1936. The theoretical 
curves of the form tV(s + 1) based on the s = 2.6 power law dependence of 
the discharge frequency distribution on discharge at this station shown in 
Fig. B. The good agreement between the theoretical curves and data shows that 
the flood recessions follow an inverse power law dependence on time which in 
turn is reflected in the power law dependence of the discharge frequency distri- 
bution. 
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Jan. Feb. Mar. Apr. May 


1955 

Figure 10. Power law fits to daily discharge data (circles) of 3 
flood recessions on the Hanochitto Rjver (7/2910) in 1955. The 

floods decay according to (t s+ ^ with s = —1.9 frccti Figure 8) . 
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Figure U. Power law fits of the form t to daily discharge 
ijg- hrt (circles) for 3 flood, decays on the Bad River (4/0270) in 
1952. Individual floods at this station indeed follow a power 
law decay, as implied by the long term discharge frequency dis- 
tribution shown in Figure 8. 
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range, from nearly 10 c.f.s. to below 10 c.f.s. 

To firmly establish the validity of our suggestion that 

the slope s of long term discharge frequency distributions 

reflects the short term decay of individual river floods 

1 

s+1 

according to t we have used the power law decays with the 
appropriate value of s, to fit ten sample flood recessions 
at each of the gauging stations of Figure 8- The results are 
summarized in Table 2 which lists, for each river in Figure 8, 
the inclusive dates of the sample flood recessions, the root 
mean square (rms) errors in the power law fits of the recessions, 
and the average rms error of fit for the 10 cases- The flood 
recessions listed in Table 2 were selected by scanning the 
available hydrologic data in the Water Supply Papers to find 
the longest uninterrupted flood decays for each river- Each 
recession was terminated when the flood had decayed to the 
stochastic level of discharge as indicated by the data. The 
rms errors of the power law fits of the 90 flood recessions 
listed in Table 2 are remarkably small, averaging only about 
5 %, a value comparable to the uncertainties in the measured 
daily discharges themselves. The smallness of these rms 
errors clearly establishes that floods decay with an inverse 
power law time dependence. 

In conclusion, we have found that floods decay with an 
inverse power law dependence on time. The ejqsonent of this 
dependence varies from river to river and even from station to 
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TABLE 2 - 


POWER LAW FITS TO RIVER FLOOD DECAYS 
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RIVER 

INCLUSIVE DATES OF FLOOD DECAYS 



H11IIIHI 

MANISTEE RIVER 

4/11-4/14/59, ,726-7/30/52, 11/21-11/24/58 

1 .7%, 7.1;', 

1 .6% 


4/1260 

4/7-4/10/56, 10/8-10/11/54, 4/11-4/14/58 

1.7%, 5.2%, 

1 .4% 

2 . 6 % 

s * - 5.0 

6/5-6/8/54, 7/14-7/18/57, 4/22-4/25/60 

0.4%, 1.1%, 

2. * 


Q«t-* 25 

5/30-6/2/53 

3.3% 



THUNDER BAY RIVER 

5/23-5/30/59, 4/7-4/1.3/59, 4/9-4/26/56 

1.3%, 3.3%, 

5.5% 


NEAR HILLMAN, MICHIGAN 

4/27-5/9/57, 3/25-4/3/53, 3/23-3/30/48 

8.7%, 2.4%, 

3 . 8 % 

4.7% 

4/1325 

10/18-10/26/54, 4/13-4/22/47, 4/9-4/15/54 

3.8%, 8.6%, 

4 . 0 % 

s = - 4.0 
Q - t-°-« 

3/15-3/21/46 

6 . 0% 



PINE RIVER 

5/21-5/30/60,. 7/10-7/20/57, 10/18-10/26/54 

6.7%, 4.3%, 

5 . 6 % 


NEAR LEROY, MICHIGAN 

6/27-7/6/54, 5/7-5/14/54, 3/24-4/2/53 

5.4%, 1.7%, 

3 . 9 % 

4.5% 

4/1250 

5/4-5/16/53, 4/8-4/21/58, 4/9-4/15/59 

1.9%, 4.7%, 

3 . 9 % 


s. = - 3 .0 
' Q « If 0,5 

11/19-11/25/58 

6.7% 



SACRAMENTO, NEAR RED BLUFF 

1/17-2/10/36, 2/24-3/23/36, 3/31-4/15/05 

•3.6%, 2.9%, 

3 . 0 % 


11/3780 

2/6-2/16/07, 4/9-5/9/07, 2/27-3/13/10 

5.2%, 8.0%, 

3 . 2 % 

4.7% 

s = - 2.6 

2/15-2/28/11, 3/22-4/1/32, 4/18-4/30/20 

6.6%, 6.2%, 

2 . 7 % 


\ Q- I' 0 - 625 

2/21-3/8/26' 

5 .5% 



WEST BRANCH OF MANISTIQUE RIVER 

4/14-4/25/53, 4/16-4/27/51, 5/C-5/21/47 

5.8%, 4.5%, 

4.1% 


NEAR MANISTIQUE, MICHIGAN 

5/5-5/14/40, 4/21-5/13/42, 4/13-4/23/56 

2.7%, 3.3%, 

5 . 4 % 

5.1% 

4/0560 

■ 1 

3/22-4/15/46, 4/30-5/20/54, 4/23-5/10/55 

6.1%, 6.9%, 

8 . 5 % 


j s = - 2.3 

! q - t -0.769 

4/26-5/9/52 

i 

3.5% 
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TABLE 2. (continued) 



RIVER 

INCLUSIVE DATES OF FLOOD DECAYS 

RMS ERROR CF FIT 

% Hm . i r 1 1 

Z D Z ^ Z 

HOMOCHITTO RIVER 
7/2910 
s = - 1 .9 
Q-f 1 - 11 

2/16-2/29/56, 5/20-6/4/53, 4/14-4/29/55 
2/7-2/14/55, 3/28-4/5/42, 2/27-3/9/53 
2/23-3/9/55, 4/18-4/27/57, 4/10-4/20/51 
4/24-5/4/51 

5.9%, 7.2%, 4.2% 
5/6%, 3.3%, 1.5% 
5.5%, 6.7%, 3.8% 
7.4% . 

5.1% 

BAD RIVER 
4/0270 
s - - 1.8 
Q - t" 1 -25 

5/15-5/28/56, 6/2-6/13/59, 10/31-11/8/55 
7/6-7/15/53, 6/4-6/13/51, 7/6-7/15/51 
5/11-5/23/50, 5/14-5/23/52, 4/23-5/5/52 
5/7-5/18/49 

4.7%, 9.9%, 6.7% 
4.4%, 5.7%, 6.7% 
3.2%, 4.0%, 4.3% 
4.2% 

5.4% 

RED RIVER 
7/3370 
s = - 1.7 
Q « t' 1 

7/19-7/26/48, 2/21-3/2/56, 2/24-3/7/51 
3/1-3/9/39, 2/7-2/15/41, 10/3-10/10/36 
4/4-4/1 4/45 , 2/26-3/8/38, 4/28-5/8/52 
.11/7-11/17/41 

5.7%, 8.0%, 8.9% 
7.7%, 8.9%, 8.1% 
■3.4%, 7.2%, 4.8% 
8.2% 

7.1% 

JUMP RIVER NEAR 
SHELDON, WISCONSIN 

5/3620 

s - - 1.5 
-2 

Q « t 

7/16-8/5/58, 5/7-5/18/58, 5/12-5/21/50 
6/11-6/21/40, 11/6-11/17/34, 11/24-12/15/34 
4/25-5/10/20, 11/13-11/24/19, 12/24-1/11/19 
11/18-11/27/26 

10.5%, 5.4%, 7-9% 
‘ 5.8%, 3. 3%, 10.5% 
8.6%, 4.8%, 8.1% 
8.3% 

7.3% 

i 
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station along the same river. Nonetheless, despite the complex 
interactions of the large number of factors which undoubtedly 
affect the flow at any point on a stream, the resultant time 
dependence of flood decays can be described by a single para- 
meter, which can be uniquely determined from long term records 
of the discharge. This power law time dependence makes possible 
the forecasting of river discharge with an uncertainty of about 
5 % for as long as a month following the flood peak. Finally 
we should note that an inverse power law dependence of the flow 
rate on time is characteristic of diffusive and random walk 
processes, suggesting a direction for future hydrologic modeling 
of the flood recession. 

2.5 Relationship Between Meander Spectra and Discharge Spectra 
The primary objective of this investigation was to 
try to determine whether significant information on stream flow 
rates could be obtained from aerial and satellite imagery of 
river meander patterns based on a possible correlation between 
the meander and discharge spectra of rivers. Such a correla- 
tion could provide the basis for a simple and inexpensive 
technique for remote sensing of the water resources of large 
geographical areas, eliminating the need for much hydrologic 
recording. 

The test of the existence of such a correlation can be 
seen in Figure 12 where we plot the median discharge Q 
determined from the discharge spectra versus the characteristic 





wave length X associated with the break in the meander 
power spectra for the rivers studied. As can be seen there 
is only a crude correlation between these two quantities, 
confirming that additional factors (see e.g. Schumm, 1968) 
also govern this relationship. This prevents a simple 
inversion of the meander power spectrum of a stream to obtain 
its discharge spectrum. The general trend of the correlation 
between median discharge and the characteristic wavelength 
is nonetheless qualitatively consistent with earlier studies 
using less quantitative methods of analysis as can been seen 
by comparison of the relationship determined by Carlston (1965). 

3. Summary 

Although NASA no longer supports this research we 
are continuing our efforts to understand the individual dis- 
charge and meander spectra and their correlations . We 
have been impressed with the uniformity in the slopes of 
meander power spectra over the small wavelength portions of 
the spectra. As X (wavelength) -»0, we have found that the 
meander power spectral density estimates behave as X a , where 
a has a value between 3 and 4 for all the rivers we have 
investigated. The value of a appears to be independent of 
discharge and other variables such as local geology, etc. 

The slope of meander power spectra for short wavelengths 
must, in some way, be characteristic of the fundamental 


meander forming processes. 

In an attempt to understand the significance of the 
essentially constant a, we are computing the power spectra 
of artificially generated meanders to see if the slopes of 
the spectra of these mathematical meanders bear any relation- 
ship to the values of a determined for rivers. In particular, 
we are testing the validity of the theory of meander formation 
proposed by Langbein and Leopold (1966). The theory asserts 
that meander formation is a random or stochastic phenomenon 
in which the meander planform is a most probable path defined 
by a random-walk model. The probability p that a river 
deviates by an angle d9 from its previous direction in a 
distance ds along the river is described by a distribution 
which is assumed to be Gaussian 


dp 


l. dglds v 
2 ' a ; 




where a is the standard deviation and c is such that ,fdp = 1. 
The solution for the most probable path, or meander form, 
between given points is the elliptic integral 

If* d0 

s 0 o . • 

«/2 (cosm-cosQ) 

where s is the arc length along the river whose local direction 
9 is measured from the line joining the end points and uu is the 
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maximum value taken by 0. 

We have computed s(9) for various values of ou and have 
numerically inverted the function to obtain 9(s). 0(s) is 

exactly what we measure in computing the meander spectra of 
rivers. We are presently computing the power spectra of 9(s) 
for various tu. Should the mathematical 6(s) yield spectra 
whose slopes are in agreement with the river value, then we 
have found an empirical verification of the minimum variance 
theory of river formation* If not, then one must question 
the relevance of the model of Langbein and Leopold to actual 
rivermeander formation. 
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APPENDIX I 


Individual River Discharge and Power Spectra 
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1 / 6109 

P W III E flE? p kS : 6o 5 IVhf?3 lH DISK KRUD5K PI50as ' f,flJ - ssan * L 



MANISTIQUE RIVER 

STATION 4/0565 MANISTIQUE. MICH. 1938-1969 
TP 4552. PIL 2. RECD5 3294-3662 PLOT 17 SET 1 


t • i i I i i i i I i ill i > i i t i i i I i i i i I 




PROBABILITY OF FLOW RATE 


MANI5TIQUE RIVER 

STATION 4/0550 BLRNEY, MICH. 1938-1969 
TP 4552. FIL 2, BECDS 2706-3071 PLOT 15 SET 1 



PROBABILITY RRTE OF FLOW LESS THAN 



h- 

for 

B 





MRNISTIQUE RIVER 

STATION 4/0495 GERMFASK, MICH. 1938-1969 
TP 4552, FIL 2, RECDS 2340-2705 PLOT 14 SET 1 






BfiBIL JTY RATE OF FLOW LESS THAN 



POWER SPECTRRL DENSITY (DEC 2 IDOO FT. 


RIVcR musk 

S9H ORTR FT5 / 90 SPEC E5T5 
1.2 X1G S FT PERCH 



1 L. 


-Li. LUJ L 


Jill 1 . 


IQ 


1 


10 


WAVE LENGTH ( IQ FT) 


1 / 3 S 79 


r hl wmi 1 B,Uh87i 1 ?« D H§K H 8aB35K F Hi3ls^a?J^aSSSe 


-( 


- 8 NG 2 




MUSKEGON RIVER 

STATION 4/1215 EVART. MICH. 1933/10-1970/9 
TP 4552, FIL 2, RECDS 3663-4106 PLOT 18 SET 1 



PROBABILITY RATE OF FLOW LESS THRN 




POWER SPECTRRL DENSITY (DEO 2 1000 FT. 


FLHM HIVth 

1136 CRTR fT3 / 90 SPEC E5T5 
1.3XICI 5 FT REACH 



i — i 1 Lul-L-i .1 i i i Li ».lx a 


WAVE LENGTH t 10 FT3 


1 / >1547 


FS BlH. S3 IN 
03:43=10.1 I 



FNG2 


PRGSRBILITY GF FLOW RR7E 


FLHHBERU RIVER RT BRBBS ISLRND 

STATION 5/3585 WINTER, WIS. 1960/10-1969/9 
TP 4552, FIL 2. RECDS 7269-7376 PLOT 31 SET 1 



PROBABILITY ARTE OF FLOW LESS THRN 



ruNcn orcLirihL ucnoi i i i ucu 


QnTQNHGQN RIVER, MICH 

1233 ORTH FT5 / 130 SPEC ESTS 
7.5 XIO^FT HERCH 



WAVE LENGTH (10 FT) 


1 / 2600 


P W BFsBW f I 



0NT0NRG0N RIVER 

STATION 4/0400 ROCKLAND. MICH. 1 942/ 1 0- 1 969/9 
TAPE 001621 RECORDS 2149-2472 PLOT 44 SET 1 








POWER SPECTRRL DENSITY (DEC 2 l OOQ FT. 


CHIP RIVER 'i' 

719 ORTR PT3 / 7S SPEC E3T3 
».7X10 4 FT RERCH 



WAVENUMBER tfCPlD^FT) 



WAVE LENGTH t 10 FT) 


1 J 5Q00 


r w awra* 3.I&P, %m , w 7zs»S6asR 


« Q>Jrr 


POWER SPEC7RRL DENSITY (DEG 2 1000 FT. 


CHIP RIVER 

743 OfiTf) PT3 / 75 3f EC E5T5 
9. 7 X ID 4 FT REACH 



HRVFNUMBEP (PER 10 5 FTJ 


10 4° 

HRVE LENGTH ( 10 FT) 

5000 / 6QCO 

r W ZS^UC?^! 1 SMfoff, °0?3>« ^SffoSK Vl303$®j5ffiA^ ONE* 



POWER SPECTRRL DENSITY (DEG 2 1000 FT. 


icrY 


10 


10' 


10 


1 0 


CHIP RIVER 


791 ORTR PT3 / 95 3PEC E3T3 


i.SXlD^FT RERfH 


V 


QjO 

© 

a 


© 

ft 


Oft 


ft 


ft) o 


a 

A 
*4. ® 




HH 16000 FT 


j 1 — i — i — lxxjlI ■ ■ ■ » < i . 1 1 


10 


-1 


10 


WAVENUMBER (PEP 10 FT) 


-J-L-i 1 * 1 LuL-1__J — I I I Ll4 

10 1 . 10° 


HAVE LENGTH t 10 FT) 


1 / 5331 







POWER SPECTRAL DENSITY (DEC 2 1000 FT. 


RIVER MNSM 

l 1030HTH PT5 / 80 SPEC E5TS 
S.3 XIO^FT REACH 



WAVENUMBER (PER10 3 FT) 



WAVE LENGTH (10 FT) 


3 / 2209 

PLI IfjCE.PSPLTl F2BLH-a IN 0I5X KfUJ03K Pr509S.HflJ.lM5Jl.L 


POWER SPECTRAL DENSITY (DEC 2 1000 FT. 





POWER SPECTRAL DENSITY (DEC 2 1000 FT. 




rivefTmnsw 

525 ORTfi PT5 / 60 SPEC E5T5 
5.6X1 oVt HERCH 


CM 

° -a 

£io 3 


UJ p 

Q io d 


o 

UJ 

CL. 

tn t 

„-io L 


- O 0 
i LJ 



4000 FT 


WAVENUMBER (PER 10 FT) 

— i In i-i-i. ,i — l i Luux 

10 1 ua° 

WAVE LENGTH 1 1CTFT) 


3400 / 5500 

P W 2§§ E aE? f io 3 fjO-lVhfS IN DI5X KflUD5K PI509S.HflJ.«?l3N.L 



SPECTRRL DENSITY (DEC 2 1000 FT. 








POWER SPECTRAL DENSITY (DEC 2 1000 FT. 1 







POWER SPECTF1HL DENSITY (OEC 2 100D FT. 




- 80 - 


MRNI5TEE RIVER, MICH 

25430STR PT5 / 200 SPEC E5T5 
4.4 X10 S FT RERCH 



10 10 3 

HAVE LENGTH t 10 FT) 


1 / 1S25S 

r hl JMW 5iMJ M BWWW*.«.« 


PR09RBILITY OF FLOW RRTE 



PROBABILITY RRTE OF FLOW LEFS THRN 





MANISTEE RIVER 

STATION 4/1240 SHERMAN, MICH. 1933/10-1970/9 
TP 4552, FIL 2, RECOS 4610-5053 PLOT 20 SET 1 


I * * i I i i i i i i i__ i i i i i i I i i » i i a * i i 





MANISTEE RIVER 

STATION 4/1235 GOATLING. MICH. 1942/10-1970/9 


TP 4552. FIL 2, HECDS 4107-4442 PLOT 19 SET 1 


* I I I t ) I » I I I I I » I « I I I • I I I I I I I 






POWER SPECTRRL DENSITY (DEC 2 1000 FT. 


FEATHER RIVER 

657 DATA FT3 / 100 3FEC E3T3 
1.2 X10 5 FT REACH 



WAVENUMBER (PER 10 5 FT) 


WAVE LENGTH C 10 FT) 


1 / 1??7 

r w as'rfrH 1 


POWER SPECTRAL OENSITY (DEG 2 1000 FT.) 


- 85 - 

H0M0CHITT0 RIVER, MISS 

2195 DATA PT5 / 200 SPEC C5T5 
2.6X10 S FT REACH 



WAVENUMBER (PER10 3 FT) 


Jm.u,j — i — i jjULLXj-a i i 

2 1 
10 10 , 

WAVE LENGTH (10 FT) 


1 / bibs 

p w bam* w 



RED RIVER ; NORTH BASIN 

STATION 540:; FARGO, N. -DAKOTA 1901-1965 
TP 4552, FIL 1 , REDDS 326-1145 j PLOT 4 SET 1; 


FLOW RATE (CU. FT. /SEC J 


- PBC33R8I! ITT RATE OF FLOW LESS THAN 



RED RIVER NORTH BASIN 

STATION 5140: FARGO, N. DAKOTA 1201-1931 
TP 4552, FIL 1, RECDS 326-690 PlOT 5 SET 1 , 

. - » ' i , f . \ 










[»!•] 



-RED RIVER; NORTH BASIN 

5'mriON 5140: FARGO, N. DRKOTR 1931-1042 
TP 4552, FIL 1, REEDS 691-816 PLBTG SET 1 - 








ruwtl -1 srtLIhHL U t n □ i IT (UtU c 1UUU M 


RED RIVER 

1297 ORTH PTS / 100 SPEC E5TS 
3.3 XlO S FT REACH 



WAVENUMBER (PER 10 5 FT) 



WAVE LENGTH t IQ ET) 


1 / S1ST 

rLT Tpr? r?rui ? * blh. 4s 

7H cas OCT CC 05tM3ilD.l 


C in 017* H 
in 0I3K KPUO 




POWER SPECTRAL DENSITY ( DEG 2 1000 FT. 


ntu tiivcn 

fMZ ORTR rT3 / 60 SPEC E3T5 
1.3 X10 5 FT RERCH 



WAVENUMBER (REF 10 ? FT) 



WAVE LENGTH (10 FT) 



1 / 3S7Q 


r W m'tPAV 8„KiE, \h VtmNmWXMlIU&M 



POWER SPECTRAL DENSITY (DEC 2 1000 FT.) 


HtL) HiVtH 

611 DRTfl FT5 / 50 SPEC E5T5 
1.3 X10 S FT RERCH 



HRVENUM8EP (PEPIo’ft) 


1X4 4-.1.. 4 1 1x4-4 -1 i A. 4 4 Li 



WAVE LENGTH (10 FT) 


£S7q / sisq 


rLT r?rUTi *6 Bit*. C in 0I3H HPU05X ri?0:?e.r1PJ.«a 

74 C35 OCT ZZ 03:43:10.1 IN 0I3K KAU03K ri303S. 


J03S.HflJ.PC0.Bl 


BAD RIVER 

STATION 11/0270 ODRNRH, MIS 
TP 4552, FIL 2, RECOS 426-681 






POWER SPECTRRL DENSITY (DEC 2 1000 FT. 


BAD RIVER 

2U10DBTH PTS / 150 SPEC Eb'TS 
i.3Xi0 5 FT REACH 



WAVE LENGTH (KTFT) 


1 / 4921 

r hl MFsEW HAhi.4" ° 15K WJOW.imj.BMIS.1 




BAD RIVER 

1051 ORTH FT5 / 100 SPEC E5TS 
5. 7X10*^ FT REACH 


PH 2000 FT 


WRVENUMBER (PER 1 (TFT) 


WAVE LENGTH ( 1CTFT) 


2700 / <1921 

PL7 TfiPE P7PL71 FI BLK. 3 IK DISK KflUmK PI3n9fi-ftfLI-Bmfl.L 


BHD RIVER, WISC 

1310 DrITR PT5 / 150 SPEC E5T5 
7,5X10 U FT REACH 


o o 

& 


& 

■% 


4000 FT 


WAVENUMBER (PER 10 FT) 


WAVE LENGTH (10 J FT) 


1 / 2S21 


p w m,irw uw«». M 


POWER SPECTRRL DENSITY (DEG* 1000 FT.) 


B 


1300 DflTI 
7 



HRVENUM8ER (PERIO'TT) 


10 1 10 ° 
HAVE LENGTH ( 10 FT) 


100 / 2700 

P W ll S IN 01 5K * flUD5K fJSOSS.rfBJ.BflDH.L 





PROBRBILITT OF FLOW RATE 


WHITE RIVER 

STATION 4/0275 ASHLANO, WIS. 1948-1970 
TP 4552, FIL 2, RECDS 682-937 PLOT 9 SET 1 



PROBABILITY RATE OF FLOW LESS THRN 







WHITE RIVER 

1030 DHTR PTS / 80 SPEC E5TS 
S.BXIO^FT REACH 


HH 4000 FT 


WAVENUMBER CPER10 FT) 


WAVE LENGTH 11 IT FT) 


2000 / 4173 


IN DISK KRUD3K PJ5.HAJ.MH3M.L 




WHITE RIVER 

2089 DATA PT5 / 125 SPEC E5T5 
I.IX10 S FT REACH 


$ 



8000 FT 


WAVENUMBER (PERIDOT) 

i- t _ hi i i t i i i In m i « 

10 l 10 ? 

HAVE LENGTH (1CTFT) 


i / mis 


P W ifSWii 5 & 83 Nhi!s tH DISK "s-WJ-w*"- 


BO I S BRULE RIVER 

STATION 4/0255 BRULE, WIS. 1942/10-1969/9 
TP 4552, FIL 2, RECDS 1-324 PLOT 7 SET 1 












POWER SPECTRAL DENSITY (DEC 2 1000 FT.) 



HflVENUMBER (PER 1 (TFT} 


WAVE LENGTH (10 FT) 


1 / W07 

P W liSSEFH 1 BLUhlfs IN DI5K KflUD5K 








B0ISE BRULE RIVER 

l»40»± DATA PT5 / 125 SPEC E5T5 
7.5 XlCJ^FT REACH 


Cl 


o 





l-H 14000 FT 





SPECTRAL DENSITY (DEG 2 1000 FT * ) 







PROBABILITY OP FLOW RATE 


RIFLE RIVER AT "THE RANCH" 

STATION >1/1395 LAPTON, MICH. 1950/10-1970/9 
TP 4552. FIL 2. RECOS 6308-6547 PLOT 27 SET 1 


I I I I « I 


PROBABILITY RATE OF FLOW LESS THRN 





RIFLE RIVER 

STRTI0N H/1405 SELKIRK. MICH. 1950/10-1970/9 
TP 4552. FIL 2. RECDS 6548-6787 PLOT 28 SET 1 







PROBABILITY OF FLOW RATE 


-107- 


RIFLE RIVER RT HIGHWAY 70 

STATION 4/1420 STERLING. MICH. 1936/10-1969/9 
TP 4552, Fit 2, RECDS 6788-7183 PLOT 29 SET 1 

±: 




RIVER RIFD 


55&BRTR PT5 / BO SPEC E5TS 


H.6 XICTFT REACH 


e m 


©a 



W2000 FT 


10 10 
WAVENUMBER CPEB10 3 FT) 


WOVE LENGTH (lO^FT) 


2750 / VVSS 


}|5 e „Ef L 3 I DISIt KWU ”* FI5096 - HSJ - I ’ IFOJ 


• i r Ji >T r i*ii ~ i~- is* 


SPECTRRL DENSITY (DEG 2 1000 FT*) 







ruNLn orcumML ucnoi i i ii/cu 


RIVER PI FM 

772 OflTfl fT3 / SO ?3T? 

l.sxi a 4 ft perch 



HRVENUMBER (PER 10 3 FT) 


HRVE LENGTH 110 FT) 


1 / 666 

W IK E BE? fL 3 J RW® IN ° 15K Km,D3H 


POWER SPECTRAL DENSITY (DEC 2 IOOO FT. I 




- 111 - 

RIVER RIFU 

12<42 ORTH PTS / 00 SPEC ESTS 
3. 3 X 1 o' 1 FT REACH 



WAVENUMBER (PER10 3 FT) 



HAVE LENGTH tlO'TT) 


1 / 12M3 

P W 39| e s|? p 5jJ 1 5i s $&S.s IN 015H KflUD3K r |5aa6 - |WJ *" iru - L 



112 - 






JUMP RIVER 

STATION 5 / 3B20 SHELDON. W1S. 1915-1970 
TP 4552. FIL 2, RECDS 7502-8151 PLOT 33 SET 1 



PROBABILITY RATE OF FLOW LESS THAN 


PRGBRBILITY OF FLOW RRTE 


-113- 


4 


SOUTH FORK JUMP RIVER 

STATION 5/3615 OGEMfl. HIS. 1944/5-1951/9 
TP 4552, FIL 2, BECDS 7377-7501 PLOT 32 SET 1 




°R0BRBIL1TY RRTE OF FLOW LESS THRN 




POWER SPECTRRL DENSITY f DEC 2 IDDO FT. 







